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EROSION  BY  CONCENTRATED  FLOW 

By  Joe  CrWillis,  research  hydraulic  engineer, 
Soil  and  Water  Conservation  Research  Division' 


When  water  flows  on  a  sloping  surface,  gravity  causes 
it  to  seek  the  lowest  levels  as  determined  by  the  topo- 
graphv  of  the  surface.  It  finds  its  way  into  the  crevices 
between  soil  aggregates,  into  furrows  and  rills,  row  mid- 
dies,  drainageways,  ditches,  small  streams,  and  rivers. 
Each  step  mav  be  looked  upon  as  progressive  concentra- 
tion of  the  flow. 

As  the  flow  is  concentrated,  its  potential  for  erosion  is 
also  concentrated.  If  the  flow  is  diffused  by  a  flattening 
of  the  slope,  its  erosive  potential  may  be  lessened,  and 
deposition  may  even  occur. 

This  report  pertains  to  the  ability  of  flowing  water  to 
entrain  and  transport  detached  material.  In  other  words. 


such  material  as  soil  particles  or  aggregates  has  already 
been  detached  and  is  available  for  transport  by  the  flow. 
Thus  only  the  transport  phase  of  the  erosion  model  as 
proposed  by  Meyer  and  Wischmeier^  is  considered.  This 
phase  was  investigated  in  a  laboratory  study  of  the  ero- 
sion of  sand  bed  material.  Laboratory  flume  tests  were 
designed  to  investigate  the  relationships  between  the  ero- 
sion rates  and  flow  variables  that  occurred  as  the  sand 
bed  was  degraded  by  set  flow  rates  and  to  define  the 
limits  of  erosion  at  the  points  of  incipient  motion  of  the 
bed  material. 


Experimental  Procedure 


Data  for  this  investigation  were  obtained  from  a  labo- 
ratory channel  16  feet  long  by  6  inches  wide  (fig.  1).  An 
adjustable  weir  at  the  downstream  end  provided  an  ele- 


In  cooperation  with  the  Mississippi  Agricultural  Experiment 
■Station  and  the  University  of  Mississippi. 


vation  control.  A  g-inch  slot  under  the  weir  passed  a 
small  part  of  the  flow  along  with  the  bed  material  and 
caused   the  sand   bed  at   the  flume  exit  to  assume  the 


Meyer,  L.  D.,  and  Wischmeier,  W.  H.  Mathematical 
simulation  of  tJie  process  of  soil  erosion  by  water.  Amer.  Soc. 
Agr.  Engin.  Trans.  12  (6):  7.54-758.    1969. 


Figure  1. -Erosion  test  channel  with  flow  near  equilibrium  for  initial  sediment  feed  of  sand  2  with  q=0.198  c.f.s.  per  foot. 


angle  of  repose  of  the  material.  This  reduced  the  effec- 
tive length  of  the  erodible  sand  bed  to  about  14.5  feet. 

The  useable  discharge  range  was  from  0.09  to  0.30 
c.f.s.  per  foot  width  of  flume.  A  vibratory  sand  feeder 
established  an  initial  capacity  transport  concentration 
greater  than  20.000  p. p.m.  A  tank  at  the  end  of  the 
channel  served  lo  trap  the  sediment  eroded  from  the 
sand  bed  and  provided  a  reservoir  for  the  pumping 
system. 

Three  sand  bed  materials  having  median  diameters  of 
0.61  (No.  1),  0.;{1  (No.  2),  and  0.40  (No.  3)  mm.  were 
used  in  this  investigation.  Figure  2  shows  the  size  distri- 
bution  of  the  three  materials.  Sands  1  and  2  were  ob- 
tained from  .sand  deposits  and  were  rather  angular.  Sand 
3  was  from  the  bed  of  Cuffawa  Creek  near  Holl) 
Springs,  Miss.,  and  its  particles  were  more  nearly  spheri- 
cal. 
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Figure  2.— Sieve  analysis  of  bed  materials. 

Sands  were  chosen  lor  this  investigation  instead  oi  silts 
or  clays  because  silt-  and  clay-sized  particles  probably 
are  seldom  detached  in  sufficient  quantities  to  limit  their 
transport  by  flow.  For  a  natural  soil,  sands  and  aggre- 
gates, whose  fall  velocities  are  adequate  to  limit  trans- 
port, comprise  a  certain  fraction  of  the  soil  mass,  and  a 
limitation  on  the  transport  of  this  material  limits  the 
erosion  of  the  soil  mass  in  general. 

Other  instrumentation  consisted  of  a  point  gage  to 
measure  elevations  along  the  bed  and  water  surface  pro- 
files, a  sonar  device  to  record  the  bed  profiles  for  rela- 
tively slow  flows  with  irregular  beds,  a  calibrated  elbow 
meter  to  determine  the  flow  discharge,  and  a  float  and 
siphon  arrangement  in  the  reservoir  to  provide  tempera- 
ture control  by  interchanging  cool  water  for  warm 
water. 


At  the  beginning  of  a  typical  erosion  test,  initial  equi- 
librium bed  and  water  surface  slopes  Were  established  by 
operating  the  flume  at  the  selected  flow  discharge  and 
sand  teed  rate.  The  erosion  processes  were  then  initiated 
by  stopping  the  sand  feed.  The  sand  bed  and  water  sur- 
face profiles  were  then  described  by  point  elevation 
readings  or  sonar  traverses  or  both  as  the  initial  sloping 
surface  was  eroded  and  degraded  by  the  constant  dis- 
charge. The  tests  were  repealed  for  I  1  to  13  discharges 
for  each  of  throe  bed  materials. 

At  the  end  of  a  typical  8-hour  erosion  test  the  bed  was 
still  being  eroded  at  a  low  rate  (about  0.003  foot  per 
hour).  Even  after  the  longest  test  used  (290  hours),  sedi- 
ment motion  was  still  observed  over  all  points  on  the 
bed.  \et  the  erosion  processes  were  approaching  some 
asymptote  or  final  coiulilion  that  would  be  reached  only 
alter  a  theoreticaliv  inlinitc  time.  Incipient  motion  tests 
were  conducted  lo  estimate  this  final  condition  for  the 
hydraulic  variables  associated  with  erosion.  The  main 
results  ol  this  phase  of  the  investigation  have  been  re- 
ported.^ 

The  procedure  for  the  final  (  onditions  tests  was  to 
allow  the  flow  and  erosion  processes  to  continue  (I )  for 
a  bed  elevation  slightly  lower  than  that  achieved  in  the 
erosion  tests,  (2)  for  the  ba.se  tailgate  setting,  (.3)  and  for 
a  discharge  giving  a  very  low  transport  rate.  The  flume 
was  allowed  to  run  under  these  conditions  for  at  least  24 
hours  to  insure  that  the  bed  forms  were  typical  of  very 
low  transport  conditions.  The  tailgate  was  then  raised 
until  incipient  motion  conditions  were  established. 

The  criterion  used  to  define  "incipient  motion"  when 
ripples  were  present  on  the  bed  was  that  motion  was 
allowed  on  not  more  than  two  nor  less  than  one  ripple 
on  a  13-foot  length  of  bed  surface.  Generally  about  15 
ripples  were  present  on  this  length  of  bed  for  the  two 
finer  bed  materials.  For  the  relatively  plane  bed  of  the 
coarser  material,  no  more  than  one  "gust"  of  sediment 
motion  over  a  2-foot  length  of  sand  bed  was  permitted 
in  a  5-second  period. 

Since  the  slopes  were  exceedingly  small  and  the  veloci- 
ties low,  subsequent  erosion  would  cause  negligible 
change  in  the  average  water  surface  elevation.  Hence,  the 
distance  that  the  water  surface  was  raised  in  establishing 
the  depth  and  velocity  of  How  for  incipient  sediment 
motion  was  taken  as  an  estimate  of  the  distance  that  the 
bed  would  be  lowered  by  erosion  after  a  very  long 
period  (theoretically  infinite  time). 


Willis,  J.  C.  Di.scussion  of  "sediment  transport  mechanics: 
initiation  of  motion,"  by  Vito  A.  Vanoni.  Amer.  Soc.  Civ. 
Engin.  Hydraul.  Div.  Jour.  93  (HYI):  101-107.   1967. 


Continuity  Equation  for  Erosion 


The   definition   of  erosion   in   terms  of  the  sediment 
transport  rate  becomes  the  continuity  equation 
dy  ^     V.G 
at  Xs 


(1) 


where  the  quantity  -^  represents  the  rate  of  change  in 

the  bed  elevation  T)  with  time  t.  or  the  volume  rate  of 
erosion  or  deposition  per  unit  area.  The  quantity  G 
represents  the  transport  rate  of  sediment  in  weight  per 
unit  time  and  unit  channel  width,  and  7^  is  the  bulk 
density  of  the  sediment,  the  quantity  V.G,  or  the  diverg- 
ence of  G,  represents  the  difference  between  the  inflow 
and  outflow  rates  of  sediment  for  an  increment  of  area. 
This  expression  is  generally  reduced  to  the  one-dimen- 
sional case  as^,  where  X  is  distance  along  a  rill,  ditch, 

dX 
or  river. 

The  continuity  equation  is  applicable  to  both  erosion 
and  deposition.  If  the  transport  rate  increases  down  a  rill 
without  a  change  in  the  flow  discharge,  then  the  rill  is 
being  eroded.  If  the  transport  rate  decreases,  then  depo- 
sition must  occur.  This  is  the  same  concept  that  Meyer 
and  \\  ischmeier"*  used  in  routing  the  eroded  sediment 
down  a  slope.   The   continuit>    relationship   is  straight- 


See  footnote  2. 
.4 


forward  ;/ we  can  predict  the  transport  rale  and  how  il 
changes  with  distance  downslope. 

Unfortunately  methods  of  predicting  the  transport 
rate  for  a  flow  arc  littic  more  than  educated  guesses. 
Several  equations  ha\e  been  proposed  for  relating  trans- 
port rate  to  flow  variables,  but  the  estimates  that  they 
yield  for  a  particular  condition  may  differ  by  more  than 
500  percent."* 

Of  course,  every  set  of  data  contains  some  error  due  to 
imprecision  of  measurement,  but  such  statistical  scatter 
could  hardly  account  for  the  observed  variations.  A  dis- 
tinction between  the  sediment  load  and  the  transport 
capacity  mav  account  for  the  extreme  variability  of 
transport  data.  Transport  capacity  may  be  defined  as  the 
rate  at  which  a  given  flow  can  transport  .sediment  of  a 
particular  particle  size  through  an  infinite  uniform  chan- 
nel reach  having  a  bed  of  the  same  material.  Hence,  the 
transport  capacity  of  a  given  set  of  hydraulic  variables 
denotes  the  capacity  of  the  .same  set  of  variables  and 
uniform  flow.  The  trapsport  capacity  then  becomes  a 
unique  function  of  the  flow  variables  and  the  sediment 
properties.  But  for  nonuniform  flow,  the  actual  sedi- 
ment load  may  differ  appreciably  from  the  transport 
capacity. 

The  actual  transport  rate  may  be  less  than  capacity 
and  this  is  illustrated  bv  figure  3,  which  shows  the  bed 
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Figure  3. -Typical  bed  profiles  for  sand  1  with  q=0.126  c.f.s.  per  fool. 
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profiles  at  three  times  during  an  8-hour  erosion  test.  An 
obvious  conclusion  that  can  be  drawn  from  this  figure  is 
that  the  erosion  rate  was  highest  near  the  upstream  end 
of  the  test  channel,  hi  fact,  0.5  hour  after  erosion  had 
begun  in  the  flume,  no  erosion  had  occurred  at  the 
downstream  end  of  the  sand  bed.  Until  this  time  the 
flow  had  apprently  entrained  a  capacity  load  by  the  time 
it  reached  the  end  of  the  flume.  At  upstream  sections 
the  load  was  less  than  capacity,  even  though  the  material 
was  sand  and  hence  was  detached  and  available  for  trans- 
port. Thus,  for  a  given  set  of  flow  conditions,  there  is  a 
maximum  rate  at  which  flow  can  entrain  a  given  sedi- 
ment. 

A  load  in  excess  of  capacity  can  be  transported  by 
steady,  nonuniform  flow.  This  is  illustrated  in  figure  4, 
which  shows  concentration  profiles  taken  30,  58,  and  82 
feet  from  the  entrance  to  a  recirculating  flume  used  in 
another  study. ^  Here  the  quantity  y  is  the  distance  from 
the  sand   bed  and   y,   is   the   flow   depth.   As   the  flow 


^Willis,  J.  C.  An  error  function  description  of  the  vertical 
suspended  sediment  distribution.  Water  Resources  Res.  5  (6): 
1322-1329.  1969. 


proceeded  down  channel  the  concentration  distribution 
became  more  nearly  in  equilibrium  with  the  flow,  and 
flow  variables  (velocity  and  depth)  changed,  causing  a 
higher  transport  capacity.  At  the  upstream  stations 
where  the  depths  were  relatively  deep  and  velocities  low, 
a  higher-than-capacity  load  was  maintained  in  the  upper 
part  of  the  flow  by  virtue  of  the  momentum  and  low  fall 
velocity  of  the  sediment  particles,  which  provide  an 
upper  limit  on  the  rate  that  sediment  can  be  deposited. 

When  erosion  occurs  in  the  field,  loads  both  less  than 
and  greater  than  capacity  can  occur.  Where  net  erosion  is 
occurring,  the  sediment  load  must  be  less  than  the  trans- 
port  capacity  because  of  the  entrainment  rate  limitation. 
Where  deposition  is  observed,  a  greater-than-capacity 
load  may  exist  for  some  distance  because  of  the  fall 
velocity  limitation  on  the  deposition  rate. 

The  continuity  equation  is  directly  applicable  to  ero- 
sion processes  once  a  capacity  load  relationship  and  an 
entrainment  function  that  describes  load  variations  are 
developed.  Until  reliable  load  and  entrainment  functions 
are  available,  the  researcher  must  resort  to  experimenta- 
tion to  generate  erosion  predictions.  This  is  essentially 
the  approach  followed  in  this  investigation. 
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4 


Sediment  Entrainment 


The  short  test  channel,  bed  irregularities,  and  small 
rates  of  change  in  bed  elevation  resulted  in  data  that 
were  only  reliable  enough  to  determine  the  average  ero- 
sion rates  for  the  test  channel.  However,  figure  3  does 
permit  a  rough  interpretation  of  how  the  erosion  rates 
varied  with  distance  down  the  channel  and  suggests  the 
form  that  an  entrainment  function  for  erosion  by  con- 
centrated flow  might  assume.  Since  the  actively  eroding 
sections  of  the  bed  profiles  can  be  approximated  by 
straight  lines,  the  erosion  rate  varied  nearly  linearly 
down  the  channel.  If  the  distance  downstream  X^  where 


a  capacity  load  would  be  entrained 
linear  relationship  becomes 


introduced,  this 


9t? 
9t 


drj 


X  =  0  ^     v/ 


(2) 


At  X  =  Xp,  the  transport  rate  is  the  capacity  rate  G^,, 
and  the  average  erosion  rate  upstream  from  this  point  is 

Q 

From  equation   (2)  this  average  erosion 


rate  can  also  be  expressed  as    ,-  r— 
^  ^  9t 

substitution  equation  (2)  becomes 


X  =  0 


3t 


~~%\\    xj 


so  that  bv 


(3) 


Integrating  equation  (3)  for  the  transport  rate  at  points 
along  the  channel  gives 


^c     ^c 


H) 


(4) 


Thus  the  erosion  rate  depends  not  only  on  the  trans- 
port capacity  but  also  on  the  relative  excess  of  the 
capacity  over  the  transport  rate.  Equations  (3)  and  (4) 
also  illustrate  the  importance  of  channel  length  in  inter- 
preting erosion  data.  Longer  channels  will  yield  lower 
average  erosion  rates,  other  things  being  equal.  Only  ero- 
sion rates  based  on  the  capacity  entrainment  distance  X^ 
will  be  unique  for  the  same  flow  variables  in  channels  of 
different  length. 

In  many  erosion  problems  the  degradation  of  the  sur- 
face may  be  negligible,  that  is,  the  overall  soil  surface 
slope  does  not  change  significantly  for  the  period  of 
consideration.  For  such  cases,  G^.  and  X^.  may  be 
approximately  constant  for  given  flow  rates.  However, 
for  the  flume  erosion  tests,  the  slope  and  flow  variables 


changed  with  time  and  the  transport  capacity  as  defined 
by  these  variables  must  have  changed  also.  In  a  long 
reach,  the  initial  degradation  at  the  upstream  end  would 
impart  observable  flow  nonuniformities,  and  transport 
capacity  G^,  and  the  capacity  entrainment  distance  X^, 
would  vary  with  both  time  and  distance  along  the  reach. 
Thus  equations  2  through  4  must  be  considered  as 
approximate  relationships  based  on  the  average  G^,  and 
Xj.  quantities  and  the  average  flow  variables  for  the  short 
test  channel.  Although  these  values  could  not  be  reliably 
determined  from  the  flume  data,  they  must  be  functions 
of  the  flow  variables  and  the  sediment  and,  therefore, 
the  resulting  average  erosion  rates  must  also  depend  on 
the  flow  variables  and  sediment  as  well  as  the  length  of 
the  test  reach. 

The  flow  variables  themselves  are  generated  by  the 
flow  quantity,  the  slope  of  the  soil  surface  or  watenvay, 
and  the  roughness  of  the  soil  surface.  Flow  variables 
such  as  velocity  and  depth  could  be  specified  but  would 
be  somewhat  difficult  to  determine  and  would  be  less 
convenient  to  apply  to  actual  erosion  problems. 

Thus    the   basic   set   of  independent   variables   upon 

which  the  average  erosion  rate  ^  or  other  variables  may 

at 

be  considered  to  depend  are  unit  width  flow  discharge  q, 
slope  S,  channel  length  L,  soil  properties,  and  fluid 
properties,  for  the  assumption  that  width  effects  are 
negligible.  In  functional  form  (f)  this  relation  becomes 

977 

g^  =  f  (q,  S,  L,  7g,  D,  s.f.,  a,  7^,  u,  g)  (5) 

where  g  is  the  gravitation  constant:  the  bulk  density  y^, 
median  diameter  D,  shape  factor  s.f.,  and  geometric 
standard  deviation  o  of  the  size  distribution  are  taken  to 
characterize  the  sand  bed  material,  and  the  unit  weight 
7^^  and  kinematic  viscosity  i>  of  the  water  specify  the 
fluid  properties. 

For  the  purposes  of  this  investigation,  since  L,  7^,  7^^ . 
u.  and  g  were  constant  or  essentially  constant,  the  func- 
tional relationship  reduces  to 


3'7_  ,■ 


9t 


I  (q,  S,  D) 


(6) 


where  D  should  be  considered  as  a  code  designation  of 
the  sediment  material  encompa.ssing  the  other  sediment 
properties  as  well,  rather  than  indicative  of  only 
diameter  dependency. 


Evaluation  of  the  Variables 


The  average  erosion  rate  -r-  is  defined  as  the  average 
rate  of  lowering  of  the  channel  bed,  which  is  equivalent 
to  the  volumetric  erosion  rate  per  unit  of  bed  area;  how- 
ever, an  incremented  approximation  as   — '    resulted  in 

considerable  scatter,  since  the  ele\ation  changes  in  short 
time  intervals  were  small.  Fortunately  the  average  bed 
elevation  was  found  to  describe  hyperbolic  functions 
with  time  of  the  form 


T?-r?„  =K{t-l^) 


(7) 


where  77^  is  the  elevation  of  the  channel  bed  at  which 
incipient  motion  of  the  sediment  would  exist  and  m,  k, 
and  tjj  are  parameters  of  the  function.  Values  of  77^  were 
determined  from  the  final  conditions  tests  as  a  function 
of  flow  discharge  as  shown  in  figure  5.  Since  only  two 
parameters  can  be  determined  from  a  least  squares 
analysis  of  the  data,  a  trial-and-error  procedure  was 
necessary  to  evaluate  the  other  parameters. 
Several  curves  of  the  form 


Vi  "Ki 


(8) 


encompassing  the  shapes  of  the  8-hour  erosion  sequences 
were  established  on  tracing  paper.  Superimposing  these 
tracings  over  the  data  curves  permitted  the  parameter 
mj  to  be  estimated  as  well  as  shifts  in  77  and  t  to  define 
the  sequence  of  TJi  and  T,  values.  Initial  estimates  of  the 
erosion  rates  were  then  made  by  differentiating  equation 
(8). 


37?^  drh 

at    at, 


T, 


(9) 


Differentiating  equation  (7)  gives  a  similar  expression 

(10) 


dr]  ^     m(r?-77o) 
at  t  -  t^ 


m  (t?-77o) 
T,  -At 


where  the  quantity  t  —  t^  is  denoted  as  a  change  At  from 

Equating  these  two  expressions  for    -^  and  grouping 

terms  gives 

^  ^1  ,3211 

T,r^-Atr^^  +  m  =  0  (11a) 


^-i?o 


^-^o 


Yi  -AtX,  +m  =  0 


(lib) 


Thus  equations  (11)  define  a  straight  line  in  the  vari- 
ables Yj  and  Xj  so  that  the  slope  -At  and  intercept  m 
could  be  evaluated  by  a  least  squares  analysis  of  the 
data.  Revised  erosion  rate  estimates  were  then  computed 
for  the  new  m  and  t  —  t^  values  and  the  least  squares 
analysis  was  repeated  until  no  additional  adjustments  in 
m  or  t  —  tp  were  indicated.  Usually  no  significant  change 
was  observed  after  the  second  trial.  The  average  erosion 
rate  associated  with  each  set  of  water  surface  and  bed 
profiles  was  then  determined  from  equation  (10)  and  for 
the  final  t^  and  m  values. 

The  slope  variable  constitutes  the  only  other  variable 
associated  with  equation  (())  that  is  not  self-explanatory. 
In  a  relatively  long  channel,  the  energv  grade  line,  water 
surface,  and  bed  will  be  essentially  parallel  and  hence 
will  define  only  one  slope  value.  For  a  short  channel  like 
the  one  used  in  this  erosion  investigation,  the  bed  and 
water  surface  slopes  may  be  different.  For  nonuniform 
conditions,  the  elevation  of  the  energ)  grade  line  H  may 
be  expressed  as 


II 


-  JL 


+  Z 


(12) 


where  y^  is  the  average  flow  depth  and  Z  is  the  water 
surface  elevation.  The  energy  slope,  which  is  the  most 
consistent  slope  estimate  for  nonuniform  flow,  is  then 
defined  as  _§H  or 

ax 


c  -  _  9H 

^ "     ax 
az 


q^  9yt 


az 


gyf  ax    ax 


(13) 


The  quantity      gv  defines  the  water  surface  slope  S^, 
and  the  depth  gradient      '  ^    is  simply  the  difference  in 

ax 

the  bed  slope  Sj,  and  water  surface  slope.  Equation  (13) 
can  then  be  expressed  as 


s  =  s,,-e3(Sw-Sb) 

g>t 


(14) 


The  bed  and  water  surface  slopes  were  then  deter- 
mined by  a  least-squares  fit  of  straight  lines  to  the 
profile  data  and  the  energy  slope  was  then  calculated  by 
equation  (14). 


I.IO 


1.06 


t 


.02  - 


0.98 


0.94 


0.90 


\  0.86 


J« 


0.82 


0.78- 


0.74 


0.70 


0.66 


0.62 


\ 


Do 


o  a  □ 


\r 


o 

03 


\ 


0.61- mm. 
bed  material 


a 


°  o.  o 


0.31-mm. 
bed  material 


0.40-mm. 
bed  material 


J L 


0        0.04        0.08 


0.12         0.16        0.20 

q.c.f.s.    per  ft. 

Figure  5. -Limiting  elevations  for  degrading  beds. 


0.24        0.28        0.32 


Experimental  Results 


The  empirical  relationships  as  indicated  by  equation 
(6)  are  shown  in  figures  6,  7,  and  8  for  the  data  from 
three  of  the  erosion  tests  with  each  bed  material.  Tests 
with  an  intermediate  discharge  are  included  along  with 
two  extreme  discharge  tests  to  sene  as  examples.  A  sum- 
mary of  the  data  from  the  erosion  tests  is  included  in  the 
appendix. 

The  data  curves  for  each  discharge  define  unique  rela- 
tionships between  average  erosion  rate  and  energ)  slope 
with  only  a  small  amount  of  scatter.  Tlie  variation  of 
erosion  rate  with  slope  and  discharge  is  clearly  illustrated 
b\   these  relationships  between  the  basic  variables.  The 


effects  of  the  different  sediments  is  not  so  well  defined 
and  seems  to  be  unimportant  except  for  small  slope 
values.  The  erosion  rate  for  the  finer  material  (0.31 
mm.)  seems  to  be  lower  than  that  for  the  other  sands  for 
similar  discharges  and  erosion  rates  of  less  than  about 
0.04  foot  per  hour. 

The  apparent  low  erodibilily  of  the  O.'il-mm.  sand  is 
supported  by  the  final  elevation  results  in  figure  5.  In- 
cipient motion  conditions  were  established  at  higher  bed 
elevations  for  this  material  than  for  the  0.40-mm. 
material.  Particle  shape  may  account  for  part  of  this 
unexpected  result.  The  grains  of  the  0.31-mm.  material 
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Figure  6.— Erosion  rate  versus  slope  and  discharge  for  ().61-mm.  sand. 
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were  quite  angular.  Also,  the  high-velocity  gradients  near 
the  bed  would  result  in  finer  particles  being  less  exposed 
than  coarser  ones  to  the  forces  of  the  flow. 

Although  relationships  like  those  in  figures  6,  7,  and  8 
are  not  dimensionally  homogeneous,  they  are  probabl) 
in  the  most  convenient  form  for  application  to  the  ero- 
sion process  in  the  field.  Discharge  could  be  determined 
from  the  excess  rainfall  rate,  and  the  surface  slope 
should  closely  approximate  the  energy  slope.  Different 
curves  for  different  soils  and  different  slope  lengths 
would  be  required,  but  our  present  knowledge  of  a  soil's 
susceptibility  to  erosion  and  of  the  entrainment  process 
is  not  adequate  to  alleviate  this  requirement. 

Similitude  principles  offer  a  means  to  make  experimen- 
tal relationships  dimensionally  homogeneous  and  inde- 


pendent of  the  scale  of  the  system.  In  connection  with 
transport  capacity  investigations  at  the  U.S.  Department 
of  Agriculture  Sedimentation  Laboratory,  Oxford,  Miss., 
the  similitude  concepts  of  fluid  mechanics  have  been  suc- 
cessfully applied  in  unifying  capacity  relationships  for 
the  available  flume  transport  data.^  The  average  concen- 
tration of  sediment  at  capacity  was  found  to  depend  on 

3 

and   a  grain   diameter 


the    Froude    number 
similitude    number 


gyt 

Similitude    theory    su^ 


lests 


^Willis,  J.  C,  and  Coleman,  N.  L.  Unification  of  data  on 
sediment  transport  in  flumes  by  similitude  principles.  Water 
Resources  Res.  5  (6):  1330-1336.    1969. 
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Figure  8. -Erosion  rale  versus  slope  and  diseharge  for  0. 40-mm.  sand. 


that  tlic  Reynolds  number  —  may  also  be  im- 
portant, but  for  the  high  Reynolds  numbers  of  the  trans- 
port capaeity  tests,  the  Reynolds  number  effeet  was  neg- 
ligible. 

The  use  of  these  similitude  numbers  in  expressing  the 
results  of  the  erosion  tests  involves  the  use  of  the  depth 
yj  in  lieu  of  the  slope  in  equations  (5)  and  (6).  Although 
the  depth  may  be  somewhat  difficult  to  determine  for 
actual  erosion  problems,  the  substitution  is  legitimate 
since  the  variables  are  interdependent. 

The  concentration  of  sediment  in  the  flow  leaving  the 
eroding  channel  was  determined  by 


3t?        L 
at  ^'  q7. 


(15) 


where  values  of  100  pounds  per  cubic  foot  for  y^  and 
14.5  feel  for  L  were  used,  fiere  C  represents  the  weight 
of  sediment  in  a  unit  weight  of  water. 

Initial  efforts  to  correlate  sediment  concentration  and 
the  Froude  number  for  the  erosion  tests  revealed  dif- 
ferent curves  for  different  Reynolds  numbers.  Dual 
Reynolds  and  Froude  dependency  can  often  be  difficult 
to  handle  for  different  scaled  systems.  Fortunately  the 
curves  of  C  versus  Froude  number  were  nearly  parallel 
on  exponential  paper,  with  the  curves  for  different 
Reynolds  numbers  being  separated  throughout  the  range 
of  the  data  by  a  nearly  constant  shift  in  the  Froude 
number.  U  this  trend  should  continue  for  all  concentra- 
tions down  to  zero,  then  physical  significance  in  terms 
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of  a  critical  Fronde  number  for  sediment  motion  could 
be  assigned  to  the  separation  of  the  data  curves. 

From  the  depth  v^.  and  discharge  q  measurements  for 
incipient  motion  testi^,  critical  Froude  numbers  for  sedi- 
ment motion  were  computed  by 

Fc-qN^  (16) 

Figure  9  shows  the  results  as  critical  Froude  number 
Fj,  versus  Reynolds  number.  Again  the  0.31 -mm. 
material  can  be  seen  to  be  less  susceptible  to  erosion  for 
low  flows  than  the  0.40-mm.  material.  For  each  of  the 
bed  materials,  the  critical  Froude  number  decreases  with 
increasing  Revnolds  number. 


figures  10,  11,  and  12  show  the  ero.sion  concentration 
as  a  function  of  excess  Froude  number  over  critical.  The 
critical  Froude  number  shifts  collapse  the  curves  for 
each  bed  material  into  one.  The  three  curves  deviate 
from  each  other  for  certain  ranges  of  F  —  F^  values; 
therefore,  the  concept  of  a  critical  Froude  number  for 
sediment  motion  does  not  adequately  account  for  the 
effects  of  sediment  properties.  Although  the  0.40-mm. 
material  was  eroded  more  easily  by  the  low  Froude  num- 
bers of  the  incipient  motion  tests,  the  concentrations  of 
this  material  that  were  entrained  during  the  erosion  tests 
were  lower  than  the  concentrations  of  the  other  two 
materials  for  most  F  —  F^.  values. 
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Figure  9.— Critical  Froude  number  for  sediment  motion  from  incipient  motion  tests. 
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Although  the  graphical  relationships  of  figures  10,  11,  ditions  meeting  this  criterion.  If  a  valid  entrainment  func- 
and  12  are  expressed  in  terms  of  valid  similitude  num-  tion  were  available,  then  results  from  an  experiment  such 
hers,  the  results  cannot  be  viewed  as  completely  general.  as  this  one  foraUmitedlength  of  erodible  surface  could  be 
If  a  capacity  load  were  entrained  before  the  end  of  the  expressed  in  terms  of  other  lengths  and  the  capacity  trans- 
channel,  then  the  results  would  be  general  for  flow  con-  port  condition  as  well. 
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Figure  10.— Concentration  of  0.61-mm.  sand  eroded  from  test  channel  as  a  function  of  excess  Froude  number  over  critical. 
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Figure  11.— Concentration  of  0.31-mm.  sand  eroded  from  test  channel  as  a  function  of  excess  Froude  number  over  critical. 
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Figure  12. -Concentration  of  0.40-mm.  sand  eroded  from  test  channel  as  a  function  of  excess  Froude  number  over  critical. 
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Summary 


Erosion  tests  were  conducted  in  a  laboratory'  channel 
to  investigate  the  relationships  betsveen  the  hydraulic 
parameters  of  shallow  flows  and  the  erosion  rates  of 
sand  bed  material.  The  results  were  expressed  in  terms  of 
the  basic  controlling  variables  of  the  flow  system  as  well 
as  dimensionless  quantities.  Although  the  basic  variable 


relationships  may  be  most  convenient  for  application  to 
similar  eroding  systems,  the  data  analysis  based  on  trans- 
port ami  flow  similitude  with  addition  of  the  concept  of 
a  critical  Froude  number  for  sediment  motion  presented 
single  curves  for  each  bed  material. 


Appendix 


List  of  Symbols 


D 
F 

f 
G 

G. 


H 
L 

m 
m, 

q 

s 

Sb 

Sw 
s.f. 

T, 
t 


-Sediment   concentration,    weight  of  sediment   in 

unit  weight  of  water. 
-Median  diameter  of  sediment  particles. 
-Froude  number. 

-Critical  Froude  number  for  sediment  motion. 
-Functional  notation. 
-Sediment  transport  rate  per  unit  of  channel  width, 

pounds  per  second  per  fool. 
-Sediment  transport  capacity,  pounds  per  second 

per  fool. 

-Gravitation  acceleration,  feet  per  second  squared. 
-Elevation  of  energy  grade  line,  feet. 
-Length  of  test  reach  downstream  of  point  of  zero 

sediment  transport  rate,  feet. 
-Exponent  of  erosion  time  variation  finiction. 

hiitial  estimate  of  m  in  trial  compulations. 

Flow  discharge  per  unit  of  chainiel  width,  cubic 

feet  per  second  per  foot. 
-Slope  of  energy  grade  line. 
-Average  slope  of  channel  b(  d. 
-Average  slope  of  water  surface. 

Shape  parameter  of  sediment  grains. 
-Initial  trial  estimate  of  quantity  (t  -  t^),  hours. 
-Time    after    initial    sediment    feed    was   stopped, 

hours. 
-Time   asymptote  of  erosion    transition    function, 

hours. 


X 
X. 


X, 


y 

>t 
Z 

7s 

Tco 

A 

V 

V 
Vi 

K 


V 


—  Distance  down  channel,  feet. 

—  Distance  down  channel  required  for  given  set  of 
flow  variables  to  entrain  a  capacity  sediment  load, 
feet.  97?  1 

—  Variabh' denoting  quantity       9Ti 

T?-T7o 

bru 

—  Variable  denoting  quantity  Ti       '      . 

V-Vo 

—  Vertical  distance  from  channel  bed,  feet. 

—  Critical  flow  depth  for  sediment  motion,  feet. 
Flow  depth,  feel. 

—  Water  surface  elevation,  feet. 

Bulk    unit    weight    of   bed    material,   pounds  per 
cubic  foot. 
-Unit  weight  of  water,  pounds  per  cubic  foot. 

—  Denotes  incremented  variable. 
Bed  elevation,  feet. 

—  Average  bed  elevation  that  should  exist  in  test 
channel  after  all  erosion  stopped,  feet. 

—  Average  bed  elevation  of  test  channel,  feet. 

—  Initial  trial  estimate  of  quantity  (V  ~  Vq)^  ff^'t- 
-Coefficient  in  erosion  time  transition  function. 

Initial  estimate  of  K. 

Kinematic    viscosity    of   water,    square    foot    per 

second. 

—  Geometric  standard  deviation  of  the  particle  size 
distribution. 

—  Laplacian  operator. 
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Table  1.— Summary  of  test  data 


Q 

Temperature 

0.10  ^ 

no 

Test 

c.f.s.  per  ft. 

^F. 

m 

Ft. 

•^c 

1-  1 

0.0926 

t    1 

0.961 

0.421 

1.068 

0.327 

1-  2 

.1256 

75.5 

1.28 

.426 

1.040 

.288 

1-  3 

.1462 

76 

1.50 

.359 

1.021 

.267 

1-  5 

.1826 

74 

1.82 

.270 

.993 

.246 

1-  6 

.2016 

73.5 

2.00 

.245 

.978 

.236 

1-  7 

.2206 

70 

2.09 

.229 

.963 

.232 

1-  8 

.2276 

73 

2.24 

.233 

.959 

.225 

1-  9 

.2690 

74 

2.68 

.204 

.927 

.211 

1-10 

.3026 

73.5 

3.00 

.175 

.901 

.203 

1-11 

.2736 

72 

2.66 

.203 

.922 

.210 

1-13 

.1724 

72 

1.68 

.300 

1.001 

.254 

2-  1 

.0950 

69 

1.12 

.199 

.955 

.140 

2-  2 

.1350 

68.5 

1.25 

.184 

.914 

.138 

2-  3 

.2008 

65.5 

1.79 

.156 

.847 

.131 

2-  4 

.2600 

68 

2.40 

.146 

.786 

.126 

2-  5 

.1980 

70 

1.88 

.163 

.848 

.130 

2-  6 

.2280 

69.5 

2.15 

.161 

.818 

.128 

2-  7 

.1800 

63.5 

1.56 

.168 

.868 

.134 

2-  8 

.2130 

65.5 

1.90 

.159 

.834 

.130 

2-  9 

.2960 

65 

2.62 

.132 

.748 

.125 

2-10 

.1118 

70 

1.06 

.182 

.938 

.140 

2-11 

.2820 

70 

2.67 

.136 

.765 

.124 

3-  2 

.1074 

70 

1.02 

.186 

.925 

.138 

3-  3 

.1290 

70 

1.22 

.179 

.894 

.130 

3-  4 

.1500 

70 

1.42 

.154 

.865 

.122 

3-  5 

.1702 

70 

1.61 

.136 

.836 

.117 

3-  6 

.1940 

70 

1.84 

.131 

.802 

.112 

3-  7 

.2080 

70 

1.97 

.123 

.782 

.110 

3-  8 

.2302 

70 

2.18 

.122 

.750 

.107 

3-  9 

.2500 

70 

2.37 

.111 

.722 

.105 

3-10 

.2738 

70 

2.59 

.106 

.687 

.103 

3-11 

.2994 

70 

2.83 

.100 

.651 

.102 
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